FTUV-99-50,IFIC-99-52 



Scalar Isoscalar pion pairs in nuclei and the A(tt, 7T7r)X reaction 

M.J. Vicente Vacas and E. Oset 
Departamento de Fisica Teorica and IFIC, Centro Mixto Universidad de Valencia- CSIC, 

4-6100 Burjassot, Valencia, Spain 
(February 9, 2008) 

Abstract 

The reaction A(ir, tttt)X has been studied at low energies, paying particu- 
lar attention to the interaction of the two final pions in the scalar isoscalar 
(I=J=0) channel. We have developed a microscopic model for the pion pro- 
duction, and then implemented the two pion final state interaction by using 
the results of a non-perturbative unitary coupled-channels method based in 
the standard chiral Lagrangians. The resulting model, describes well the 
reaction on the nucleon for all different isospin channels. Finally, we have 
considered the reaction in nuclei. Our calculation takes into account Fermi 
motion, Pauli blocking, pion absorption, and also the strong modification of 

the 7T7T interaction in the nuclear medium. 
13.60.Le,13.75.Lb,13.75.Gx,21.65+f,25.80.Hp 
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I. INTRODUCTION 



The reaction A(tt, titt)X has attracted much interest lately, as a means to investigate the 
properties of correlated I = J = pion pairs, the "<r" meson, in nuclear matter fflfl]. This 
has been partly motivated by recent experimental results, which show an A dependent large 
enhancement of the cross section at low invariant masses of the dipion system The 
enhancement occurs for the A(n~ , tt + tt~)X process (I n7T = and 2 at low energies), but is 
conspicuously absent for the A(tt + , -k + -k + )X case (I n7T = 2). Thus, it happens only when 
I wn = pion pairs are produced, and could be related to the medium dependence of the nn 
interaction in the a channel. 

Our theoretical understanding of the tttt scattering in vacuum has improved considerably 
in the last few years. Close to threshold, the nir interaction is well described by Chiral 
Perturbation Theory (xPT). At higher energies, where unitarity matters and we are out of 
the %PT range of validity, some very successful non perturbative models have been developed 
|5|||. They are able to give a good description of the tctt scattering in the a channel up to 
energies above 1 GeV. 

Whereas in vacuum the low energy tttt interaction in the I = J = channel, although 
attractive, is too weak to admit a bound state, it was soon realized that the nuclear medium 
attraction of the pions could lead to the accumulation of some strength close to the two 
pion threshold, or even to the appearance of a new bound state f7|. A similar conclusion 
has been reached in Ref. f2j, where an enhancement of the the spectral function in the a 
channel, just above the 2n threshold, has been predicted as a consequence of a possible 
partial restoration of chiral symmetry in the nuclear medium. Actually, in more detailed 
calculations, based on the irn interaction model of ref. || it has been found that spontaneous 
s-wave pion pair condensation could appear at densities as low as po [|SHT0|j. Further works 
along the same line, but using the chirally improved Jiilich model for the meson-meson 
interaction, have shown that instabilities are pushed up in density when chiral constraints 



are imposed [11, 12]. Nonetheless, some accumulation of strength close to the two pion 
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threshold was found. Similar results have been obtained in ref. JT3|]. The latter work is 
based in the model for irn interaction of ref. ||, which generates the scattering amplitude 
using the chiral lagrangians, and then unitarizes it by the inverse amplitude method. 

The first experimental signals showing a small enhancement of the irn spectral function 
in the A(ir, 2ir)X reaction were found in ref. [H], although the experiment didn't measure 



below M nn = 300 MeV, and hence, missed the most interesting region. The advent of 
CHAOS improved considerably the experimental possibilities, and it has provided us with 
a wealth of good quality data for several charge channels and nuclei. The results show the 
emergence of a very clear peak in the spectral function close to threshold, when pions 
in the scalar isoscalar channels are present in the final state P,|^JT5[|. This peak, absent 



for the reaction on deuterium already appears on light nuclei, and grows as a function of 
the atomic number. In contrast, only minor changes are observed when the cross section 
of the (tt + , 2tt + ) reaction in deuterium is compared to the same process in heavier nuclei. 
Certainly, the different behaviour of the two processes cannot be caused by trivial nuclear 
effects, like Fermi motion, Pauli blocking, or pion absorption, which are practically identical 
for both cases. 

In ref. a calculation of the A(n, 2n)X reaction has been presented in which the irir 
FSI was taken into account using the approach of ref. fl!] . The emphasis being placed 



on the analysis of FSI, the elementary mechanism of pion production on a single nucleon 
(ttN — * tcttN) was described in a simplified manner, taking only the most important pieces of 
the amplitude. Also, some simplifications were done on the treatment of nuclear effects, like 
Fermi motion, or pion absorption. The results are very stimulating, and show some enhance- 
ment on the M n7T distribution close to threshold, which agrees well with the experimental 
data. 

Our aim in this paper is to do a more detailed study of the reaction, including a realistic 
ttN — > ttttN amplitude, and incorporating some nuclear effects omitted in the previous work, 
and which could affect its conclusions. 

In sec. [n] we give a brief outline of the model for tttt scattering in the scalar isoscalar 



channel both in vacuum || and in nuclear matter [pl|. Then, in sec. [Ill] we develop a model 
for the ttN — > rnrN reaction. Finally, in sec. [TV] we discuss the process in nuclei. 



II. BASIC THEORY OF THE vrvr SCATTERING IN THE SCALAR ISOSCALAR 

CHANNEL 

A. 7T7T scattering in vacuum 

The 7T7T scattering amplitude is calculated solving the coupled channels Bethe-Salpeter 
(BS) equation 

= Vvtjt + V-nlGuT^, (1) 

where the subindex ir corresponds to the state \tttc,I = > and the subindex I accounts 
for the \tttt,I = > and the \KK, I = > states. The potentials Vy are obtained from 
the lowest order chiral Lagrangians [16| . Explicit expressions can be found in Ref . []13| . The 
VQT term of eq. [I] stands for 

VQT = J 7^4 V(ft, q 2 , k)G(P, k)T(k; q[, q' 2 ), (2) 

where q, q', P, and k are the momenta of the mesons as defined in fig. [I]. A cutoff (q ma x = 1 
GeV) is used in the evaluation of the integral. Q(P,k) is the product of the two meson 
propagators in the loop, 

<3^ k )= l k 2_ m 2 +le (p_ k) 2_ m 2 + te ( 3 ) 

The particular Vij off-shell dependence leads to the simplification of the integral BS equation 
which can be transformed into the purely algebraic expression 

Tik = Vik + VijGjjTjk , (4) 

where Vu is now the on-shell part of the potentials (when the momenta are taken such that 
pf = m]), and 



33 l J (2vr) 4 k 2 - m\ + xe (P - kf - m| + «e ' 1 J 

The resulting amplitude reproduces well experimental phase shifts and inelasticities in the 
I = J = channel up to 1.2 GeV. 



B. 7T7T in nuclear matter 

In the nuclear medium, the pion propagators are strongly modified, mainly due to their 
coupling to particle-hole (ph) and A— hole (Ah) excitations. Therefore, new terms, of the 
type depicted in fig. § have to be included in the calculation of the scattering amplitude. It 
has been shown [|17| , |l"3f that the terms c),d),f),... of fig. Q cancel the off shell contribution 
of the terms a),b),e),... so that only these latter terms must be considered with the nir 
amplitudes factorized on shell. As a result the new BS equation is written as 

T n7T = K-tt + VttkGkrTktt + Kr7r^7T7r^7T7r- (6) 

As the KK state contributes very little to T nn at low energies, the vacuum values are kept 
for Gkk- Only T nn and G nn are renormalized in the medium. The function G^ is now 
given by 

r - • ( 1 1 (7) 

** ~ J (2tt) 4 k 2 — ml- U(k) (P - k) 2 - m) - U(P - k) ' [) 
where II(/c) is the pion selfenergy in the nuclear medium calculated as 

(£) 2 fcM*0 
2 g'U(k) 



n( fc ) = ,/fV r .rrL > W 



with g' the Landau-Migdal parameter (g' = 0.7), / the iViW coupling constant (/ = 1), 
and U the Lindhard function for both ph and Ah excitations ||18|| . In addition, we have also 
included the contribution of 2p2h excitations as explained in ref. fl"3|j . 

In fig. ^| we show the imaginary part of the nir scattering amplitude for several nuclear 
densities. The results show a displacement of strength towards low energies. Although 
eq. |] leads to a zero value for the amplitude close to the Adler's zero, (s = m 2 /2 with s 



the Mandelstam variable of the two pions system), the accumulation of strength close and 
below the two pion threshold is large and grows rapidly as a function of the nuclear density. 
Very similar results were found in other works using different models for irn interaction but 
imposing some minimal chiral constraints. See, for instance, fig. 12 of ref. |12| and fig. 3 of 
ref. ra. 



III. THE ttN -> vrvriV REACTION 

Our understanding of the nN — > ttttN reaction at low energies has improved considerably 
in the last few years both experimentally and theoretically ||20-p4|. The model we present 



here follows closely that of ref. p0| , where only the 7r p — > 7i + n n process was calculated. In 



this paper we will extend it to other charge channels and will include additional mechanisms, 
that have been found to be important in the calculation of some differential cross sections 

Hi- 

At low energies, the only particles that need to be considered are pions, nucleons, and the 
A and Roper resonances. The need of the latter seems surprising, because of its relatively 
high mass and the smallness of the irNN* coupling. However, the Roper resonance partly 
decays to a nucleon and a s-wave pion pair. This mechanism does not vanish at threshold, 
unlike most other terms, and it has been found to be very important for all charge channels 



that allow for an I=J=0 pion pair [25 



We have considered in our model the mechanisms represented by the diagrams of fig. |j. 
The dashed lines are pions, and the internal solid lines are all possible baryons (N, N*, A). 
Notice, that different charge channels allow or forbid some of the diagrams. Also all possible 
different time orderings, not explicitly depicted, are included in the actual calculation. The 
Feynman rules used to calculate the scattering amplitude are derived from the effective 
Lagrangians of the Appendix. It is relevant to mention here that no attempt to fit the 
ttN — > ttttN data has been made. Rather, standard values for the coupling constants, 
obtained from different experiments and/or analysis, have been used. 
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A. Isospin amplitudes 



The scattering amplitudes for the different irN — > ttttN charge channels can be written as 
linear combinations of four isospin amplitudes T21J . Here, / and I n>n are the total isospin 
of the system and the isospin of the two final pions. The expressions, for the channels we 
are interested in, are 

T^-p - vr+vr-n) = -4=T- i2 + ^T 10 + \{T n - T 31 ) , (9) 



T (n-p -> vrVn) = -i=T 32 + ^T 10 , (10) 

T ^+p ^ ^+ 7r + n ) = -4=T 32 . (11) 

Before proceeding, it is interesting to make some qualitative considerations. The symmetry 
of the pion pair wave function implies that the orbital angular momentum of the two final 
state pions will be even for the T 10 ,T 32 and odd for the 7 T7r = 1 amplitudes. Thus, at low 
enough energies, where only s-wave would dominate, both T\\ and T 31 should be negligible. 
On the other hand, the cross section for the n + p — > 7i + 7i + n channel is much smaller than for 
the other two cases. Therefore T 32 < T 10 , what implies that the amplitudes for the 7r + 7r~ 
and for the 7r°7r° case should be equal at low energies and also that they are, essentially, 
dominated by the production of scalar isoscalar pion pairs. As a consequence, the 7r°7r° cross 
section should be about one half the one for tc + tt~ , once the different thresholds are taken 
into account. Of course, the energies at which there are experimental data are not that 
close to threshold and one cannot neglect T\\ and T 31 . In our calculation, we will proceed 
as follows. First, using the Lagrangians of the appendix we calculate the amplitudes for the 



different charge channels. Then, by means of eqs. |9], [TO] and [11], we obtain T 10 , T 32 , and the 
combination T n — T 31 , and finally, we modify T 10 by including the final state interaction of 
the I = J = pion pairs . 

At low energies, due to the smallness of the I n7T = 1 contributions, we only need to 
consider the effects of the FSI on the I nn = 0,2 channels. Furthermore, as it is shown in ref. 



IJ the 7T7T interaction in for I nn = 2 is quite weak and does not change appreciably inside 
the nuclear medium. Thus only the FSI on the I n7T = channel could produce large effects 
on the scattering amplitude both in vacuum and in nuclei. To account for it, the amplitude 
Tio is modified in the following manner, 

Tio = T\^T = Tio + TioG n7r T nn (12) 

where G n7T and T n7T are the two pions propagator and the scalar isoscalar two pions amplitude 
defined in eqs. and || respectively. Using again eqs. |9|, [U] and |Tl|, where T w is now replaced 



by Tio we obtain the scattering amplitudes for the physical channels. 

As can be seen in fig. this model, without free parameters, reproduces fairly well the 
total cross sections for the three channels considered. This agreement extends to differential 
cross sections, as it will be shown for some cases in next section. The net effect of the 
inclusion of the pions FSI is a small enhancement of the total cross section for the n + ir~ and 
7r°7r° channels, and slight modifications of the differential cross sections, as it was expected, 
given that the T factor changes very smoothly, in vacuum, over the available phase space. 

IV. THE -it A -» ttttX REACTION 

There are several nuclear effects that could modify he pion production cross section, 
like Fermi motion, Pauli blocking, pion absorption and quasielastic scattering. Additionally, 
there could be new reaction mechanisms involving more than one nucleon, but the contri- 
bution of this kind of processes has been shown to be quite small ||. To account for the 
medium effects we follow ref . [B3J . Assuming only one nucleon mechanisms, the cross section 



can be written as 



a = IfrtdzFm J |* / H / - n(W + k - ,1 - (13) 

E l T | 2 7^ (g° + E (k) - - - E(q + k-qi- q 2 )) F out _ abs , (14) 

M^i) 2w(g 2 ) 

where q,k,q[, and q~2 are the initial pion, initial nucleon and two final pions momenta; 
E(k),uj(q[) and uj(q~2) are their energies. The spatial volume element is written in terms 
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of b, perpendicular, and z parallel to the beam momentum. The occupation number n(.) 
refers to the local density, it takes the value 1 when the argument is below the local Fermi 
momentum (fcp = (3n 2 p(b, z)/2)) and zero above it. Note the sum over initial and final 
spins, but not over isospin because the channels we study involve only protons in the initial 
state and neutrons in the final one. In this paper, we only consider symmetric nuclei, and 
we will use the same density for neutrons and protons. The nuclear densities are taken 
from ref. |j34| . The amplitude T is also evaluated at the local density, what fundamentally 
means that its only changes occur in the T factor, related to the titt interaction. Finally, the 
factors Fjsi, and F out _ a b s account for the flux loss due to pion absorption and scattering in 
the case of Fjsi and for pion absorption alone in the case of F out _ a bs- Both pion absorption 
and quasielastic scattering are quite strong at the energy of the incoming pion and both 
reduce the effective initial pion flux. That is clear for the absorption, but it is also true 
for the 7r-nucleus quasielastic scattering, because in these collisions the pion loses always 
some energy, what reduces enormously the possibility of a subsequent pion production. We 
implement the flux lost with the eikonal factor 

F ISI = exp ( / dz' (P abs (uj(q),p) + P qua (uj(q), p)), (15) 

J — oo 

where P a bs,Pqua are the absorption and quasielastic scattering probability per unit length, 
that we take from ref. |33] . Due to these interactions, the flux reaching the inner nucleus is 



quite small and the reaction happens at the surface, what reduces considerably the possibility 
of strong medium effects. As mentioned in ref. Jl|, the final pions have low energy and 
absorption or quasielastic scattering is a minor effect in their case. This has also been 
confirmed in the analysis of experimental data of ref. jnj] . Nonetheless, it could affect more 
those events happening at high densities in the center of the nucleus. Therefore, we include 
the factor F out _ abs that calculates the reduction of the cross section due to the absorption of 
any of the final pions. It is given by 

/OO pOO 
_ dh P a bs(v(qi),p)) exp ( / ^ dl 2 P a f>s(w(g 2 ), p)) , (16) 
-b,z J —b,z 



where dl\, dl 2 are the elements of longitude along the <fi and q 2 directions. The probability 
of quasielastic collisions is very small at these energies and we ignore it. 

V. RESULTS 

The results that will be presented in this section are all compared with CHAOS data. 
We have approximately implemented the experimental acceptance cuts in our calculation. 
Let us take the beam direction as the z— axis, and x forming with z the horizontal plane, 
then cf) = ± 7 degrees, where <fi is the pion angle with the xz plane, and 10 < 9 < 170, with 
9 the pion angle with the z direction. Also, we have taken a pion kinetic energy threshold, 
T w > llMeV. 

In fig. H we show the invariant mass distributions for the deuteron case. We obtain a 
quite good agreement in the 7r + 7r + channel in both shape and size. The agreement is not as 
satisfactory for the 7r + 7r~ case. As can be seen in fig. [5], the total cross section is slightly 
overestimated by our model, and a similar situation is found for the mass distribution. 
In order to compare easily with the experimental shape, we have renormalized our result 
reducing it by a 20% and this is what is shown in the figure. Let us remember here, 
that no parameter has been adjusted, and some of them, like those related to the Roper 
properties are very uncertain. A quite small change of a few percent in the Roper coupling 
constants would produce a fine agreement in size. It is interesting to understand the quite 
different behaviour of the two channels under consideration. As it was already stressed in 
ref. H, the tt + tt + case follows closely a pure phase space distribution. The two peaks of the 
figure respond only to the geometry of the experimental apparatus, which favours clearly 
the situations in which either the pions go together (low M nn ) or in opposite directions 
(high M n7T ). The smallness of the distribution in the tt + ii~ reaction reflects the much 
richer structure of the amplitude, and it is produced by the destructive interference of large 
pieces. Consider, as an example, the mechanisms (I) ir~p — > n* — > n(iTTc)j = j = o and (II) 
7i~p — > n* — > An — > ri7T + 7T~. The (I) amplitude has a constant sign, all over the available 
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phase space. However, process (II) has an amplitude approximately proportional to the 
scalar product of the final pions momenta. Therefore, it changes sign when passing from low 
to high invariant masses. In fact, not only is there destructive interference at low masses, 
the constructive interference increases the size of the high mass region. 

Our results for Calcium are shown in fig. [7|. Again, there is a very good agreement 
with the 7r + 7r + data. The shape of the figure is quite similar to the deuteron case. There is 
some softening, and the distribution reaches higher masses. Both features are mostly due to 
the Fermi motion of nucleons. Pion absorption is weak at low energies. Hence, it does not 
affect much the final pions, although it is partly responsible for the reduction of the spectral 
function at high masses. More important is the initial pion absorption, which apart from 
changing the total cross section, prevents the pions from reaching the nucleus core, therefore 
decreasing the effective density at which the reaction happens, and modulating all nuclear 
effects. 

The quality of the agreement on the latter channel, gives us confidence that all stan- 
dard nuclear effects are properly taken into account. In particular that we have a proper 
description of where in the nucleus the reaction takes place. 

Finally, let us discuss the ir + 7i~ data, the channel in which the nuclear dependent tttt 
interaction in the scalar isoscalar channel could be responsible of the presence of a low mass 
peak. Our results show some enhancement close to threshold, even when the pion FSI is 
calculated in vacuum. This is due to Fermi motion. In the full calculation, the inclusion 
of the medium dependence in the FSI leads to a further enhancement and a displacement 
towards low masses, although is not enough to reproduce the data. 

One first question is the consistency of our results with those of ref. which repro- 
duce well the experiment. In order to do a proper comparison, we should impose some 
approximations that were used there. In particular, we find that pion absorption forces the 
reaction to occur at the nuclear surface, at densities lower than those used in ref. [I]. If we 
impose a high fixed average density (see fig. |^), we also get a large, although insufficient 
enhancement, and too high values for the mass distribution at high masses, produced by the 
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excess of Fermi motion. If we neglect pion absorption, the results follow closely the curve 
corresponding to p — 0.5p , as expected, given the fact that for the Calcium nucleus, the 
average density is approximately half the nuclear density. 

In order to facilitate further the comparison with the results of ref. |1[ we have also 
considered a simplified model for the reaction, neglecting all terms with a A resonance. The 
two curves in fig. |9] correspond to FSI calculated in vacuum and with a density p = 0.5po- 
This simplified model has a similar structure to that of ref. Q, and although it is able to 
reproduce the peaks structure in nuclei, it does so at the price of overestimating for the 
deuteron case the threshold region by a large factor. 

VI. CONCLUSIONS 

We have studied the A(tt + , 7r + 7r )X reactions on the deuteron and on Calcium using a 
realistic model for the elementary tcN — > ttttN production, and including several nuclear 
medium effects, like Fermi motion, pion absorption, pion quasielastic scattering, and the 
medium dependent tttt interaction in the scalar isoscalar channel. 

We find a very good agreement with the experimental data for the ir + n + production, 
which gives us confidence on our treatment of the common nuclear effects. However, we are 
unable to reproduce fully the strong enhancement close to the two pion threshold found in 
the experiment for the 7r + 7r~ production. We also find that approximations previously used 
in the literature could be critical in reproducing such an enhancement. Several possibilities 
open up. It could happen that the tttt interaction in the a channel has a stronger depen- 
dence on density than that provided by existing models. This would be most interesting. 
However, we have also found that the smallness of the spectral function in the deuteron, 
close to threshold, is due to destructive interference between large pieces of the amplitude. If 
some of these pieces is substantially modified in the nuclear medium, the interference could 
disappear, and the spectral function would look more like phase space, as it is the case for 
7r + 7r + production. This would be enough to reproduce the experimental data. New dedi- 
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cated experiments, with a wider phase space, and at lower energies, where the interference 
effects are smaller are important to settle these questions. 

We also find that absorption of the incoming pion, leads the reaction to occur at the 
surface, therefore reducing the signals of any density dependent effect. Electromagnetically 
induced reactions, free from such a disadvantage are clearly called for to investigate the 
medium effects on the tttt interaction in the a channel 
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APPENDIX A: LAGRANGIANS USED IN THE irN -> vrvriV MODEL 

1. Pions and nucleons 

The scattering amplitude of the mechanisms depicted in fig. [|, considering only nucleons 
for the internal lines, can be derived from the chiral lagrangians with the inclusion of baryons 



of ref. ||T6| , |35| -p^1 . The relevant pieces can be written as 

L = L W7r + L nN , (Al) 

In this equation, L nn contains the purely mesonic interaction and L n N the meson(s)-nucleon 
terms. The mesonic part is given by 

Anr = ^ < dptfdHJ + X (U + U^) >, (A2) 

The brackets indicate the sum in flavour space. The matrices U and u are defined by 

U($) = w($) 2 = exp {iV2$/f} (A3) 

where 
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\ 



V2 W ) 



(A4) 



where r are the Pauli matrices and the w are the pion fields. At lowest order, / is equal 
to the pion decay constant, / = f n — 92AMeV. Assuming isospin symmetry, the mesonic 
mass matrix is 



(ml ^ 



X 



V 



ml j 



(A5) 



The interaction with the nucleon is described by the term 



9a 



Here, M is the nucleon mass, ^ is the nucleon isospinor 



(A6) 



( \ 

P 



(A7) 



and qa ~ 1-26 is related to the pion nucleon coupling constant (/jvjvttAtv) by the Goldberger- 
Treiman relation = |j . The covariant derivative of the nucleon field V^\l/ is given by 



(A8) 



Expanding U and w, and keeping the terms with up to four pion fields, the following set 
of Lagrangians are obtained, 



1 



j NNtv 



fNNn 



(drf$) 2 -p(d lt $) 2 + -mlP 



m, 



j NNtttt 



6/2 
1 



^7 M 75[(^0r)0 2 -(0r)(^00)]M/ 



(A9) 
(A10) 
(All) 
(A12) 
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Whereas the expression of the Lnn-k Lagrangian is quite standard, different forms can 
be found in the literature for the |A9| and [All| pieces, because different representations for 



the pion field U(<&) had been used. However, it can be shown that they produce the same 
physical amplitudes. In particular, they are equivalent to the Weinberg Lagrangians used 



in ref. |2tJ when the chiral symmetry breaking parameter £ is taken to be zero. 

The Lmn-k-k term describes the isovector part of the ttN — > ttN s-wave interaction. The 
very small isoscalar part would appear at higher orders of the chiral Lagrangian. Adding 
a phenomeno logical isoscalar part, and doing a low energy approximation, the Lagrangian 
can be recast in the typical form 



wn™ = -4?r j^0 2 ^ + ^ff(M0)*l . (A13) 
\m n mi I 



where A2 = 167r j 2 ~ 0.045. Fitting the constants to the s-wave nN scattering lengths p8 
one gets an slightly larger A 2 = 0.52, and Ai = 0.048. 
Finally, the pion and nucleon propagators are 

DM = 2 1 2 , (A14) 
q 2 — m 2 

= (A15) 

where E q = y/q 2 + M' 2 

2. A resonance 

Diagrams d) and e) of fig. |] could also have A resonances as intermediate steps. The re- 
quired AA^ and AA7r vertices are described by the following phenomenological Lagrangians 



m 



A/44(^) t Vn + h.c. , (A16) 



Laa« = —ASAii^T^A , (A17) 
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where the ^aw) are two-component spinor fields, (T*), and 5a (7a) are the spin (isospin) 
1/2 — ► 3/2 and 3/2 — > 3/2 transition operators. Definition, and some useful algebraic 
relations of this operators can be found in ref. j|39|| . 
The A propagator is given by 

G A (p) = — (A18) 

where W is the A invariant mass and the resonance width Ta is 

1 / f* \ 2 M 

r *M = *;{k) W^W-M-m.). (A19) 

with q cm the pion momentum in the resonance rest frame. The AiW coupling constant, 
obtained from ttN phase shifts or from the experimental width, using eq. |A19| , takes the 



value /* = 2.13. For the AA7r coupling we take the quark model value, f A = 4:/5fNNw [40 



3. Roper resonance 

The Roper resonance can decay into a nucleon and a pion, a A and a pion and a nucleon 
and two s-wave pion. The effective Lagrangians we use to describe these interactions are 

f - 

L n *n, = -—^ N *i^d^Ty N + h.c. , (A20) 

m n 

L N * A . = f -^^^ A S\{dS)TH N , + h.c. . (A21) 

and 

L N *Nnn = -CV N *$$if>N + h.c. . (A22) 

The coupling constants /, /at*Att and C are calculated from the A^* width and branching 
ratios. There are considerable uncertainties in the experimental information, but using the 
central values from ref. [[|l]] one gets / = 0.477, /at* Att = 2.07 and C = 2.3m~ 1 []39|| . The 
Roper propagator is given by 
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G - (P) = W-M-Uw N , iP) (A23) 
At low energy, the width is dominated by the decay into the Nir channel, and it is given by 

3 / f \ 2 M 

T N * « T Nn = — M- — |q cm | 3 Q(W — M — m n ) . (A24) 
2n \m n J W 

Whereas all other constants are obtained from different experiments, without two s- 
wave pions in the final state, and hence are independently determined, we should handle 
with more care fN*An and C. Both constants correspond to processes with two pions in 
the final state, which could be in the scalar isoscalar channel. The constants have been 
fitted to the experimental Roper partial widths using the Born approximation (without 
final state interaction of the pions), and therefore effectively incorporate already the vacuum 
renormalization due to the final state interaction of the pions. As we will include the FSI 
explicitly in our calculation, we have first to discount its effects on these two constants to 
avoid double counting. 

As explained in Sec. |111 A| our model implements FSI by multiplying the scalar isoscalar 



part of the amplitude by the factor T. This factor depends smoothly in vacuum on the 
invariant mass of the two pion system, and is practically constant over the available phase 
space region of the AT* (1440) decay. We thus take a corrected value for C 

C ^wy (A25) 

where T is the value of T at the average invariant mass of the two pions in the A/"* (1440) 
decay. The correction is not as important for the fN*-*An case. First, because only a small 
part of the partial width produces pion pairs in the scalar isoscalar channel, and second 
because at the low energies of the experiment under analysis, the contribution of this kind 
of mechanisms is small. 



17 



REFERENCES 

[1] R. Rapp et al, Phys. Rev. C59, R1237 (1999). 

[2] T. Hatsuda, T. Kunihiro and H. Shimizu, Phys. Rev. Lett. 82, 2840 (1999); T. Hatsuda, 
T. Kunihiro, H. Shimizu, |nucl-th/ 98 10022] T. Kunihiro, |hep-ph/ 9905252] . 



[3] F. Bonutti et al, Phys. Rev. Lett. 77, 603 (1996). 
[4] F. Bonutti et al, Nucl. Phys. A638, 729 (1998). 

[5] G. Janssen, B.C. Pearce, K. Holinde and J. Speth, Phys. Rev. D52, 2690 (1994); 
D. Lohse, J.W. Durso, K. Holinde and J. Speth, Nucl. Phys. A516, 513 (1990). 

[6] J. A. Oiler, E. Oset, Nucl. Phys. A620, 438 (1997); J.A. Oiler, E. Oset and J.R. Pelaez, 
Phys. Rev. Lett. 80, 2452 (1998). 

[7] P. Schuck, W. Ndrenberg and G. ChanfrayZ. Phys. A330, 119 (1988). 

[8] G. Chanfray, Z. Aouissat, P. Schuck and W. Ndrenberg, Phys. Lett. B256, 325 (1991); 

[9] V. Mull, J. Wambach and J. Speth, Phys. Lett. B286, 13 (1992). 

[10] Z. Aouissat, G. Chanfray and P. Schuck, Mod. Phys. Lett. A15, 1379 (1993). 

[11] R. Rapp, J.W. Durso and J. Wambach, Nucl. Phys. A596, 436 (1996). 

[12] Z. Aouissat, R. Rapp, G. Chanfray, P. Schuck, and J. Wambach, Nucl. Phys. A581, 
471 (1995). 

[13] H.C. Chiang, E. Oset and M.J. Vicente- Vacas, Nucl. Phys. A644, 77 (1997). 
[14] P. Camerini et al, Nucl. Phys. A552, 451 (1993). 
[15] F. Bonutti et al, Phys. Rev. C (1999) 

[16] J. Gasser and H. Leutwyler, Nucl. Phys. B 250, 465 (1985). 
[17] G. Chanfray and D. Davesne, Nucl. Phys. A 646, 125 (1999). 



[18 
[19 
[20 
[21 
[22 

[23 
[24 
[25 
[26 

[27; 

[28 
[29 
[30 
[31 
[32 
[33 



E. Oset, H. Toki, W. Weise, Phys. Rept. 83, 281 (1982). 
P. Schuck et al. , frucl-th/ 9806069 . 



E. Oset, M.J. Vicente Vacas, Nuc. Phys.A 446, 584 (1985). 

O. Jaekel, M. Dillig, C. A. Z. Vasconcelos, Nucl. Phys. A541, 673 (1992). 

V. Bernard, N. Kaiser and U.G. Meissner, Nucl. Phys. B457, 147 (1995); ibid. A 619, 
261 (1997). 

T.S. Jensen and A.F. Miranda, Phys. Rev. C55, 1039 (1997). 

A.A. Bolokhov, M.V. Polyakov and S.G. Sherman, Eur. Phys. J. Al, 317 (1997). 

V. Sossi et al, Nucl. Phys. A 548, 562 (1992). 

G. Kernel et al, Phys. Lett. B216, 244 (1989). 

C.W. Bjork et al, Phys. Rev. Lett. 44, 62 (1980). 

J. Lowe et al, Phys. Rev. C44, 956 (1991). 

A.A. Belkov et al, Sov. J. Nucl. Phys. 31, 96 (1980). 

S.A. Bunyatov et al, Sov. J. Nucl. Phys. 25, 177 (1977). 

M.E. Sevior et al, Phys. Rev. Lett. 66, 2569 (1991). 

G. Kernel et al, Z. Phys. C48, 201 (1990). 

E. Oset, M.J. Vicente Vacas, Nucl. Phys. A 454, 637 (1986). 



[34] C.W. De Jager, H. De Vries and C. De Vries, Atom. Data Nucl. Data Tabl. 14, 479 
(1974). 

[35] U.G. Meissner, Rep. Prog. Phys. 56, 903 (1993); V. Bernard, N. Kaiser and U.G. 
Meissner, Int. J. Mod. Phys. E 4, 193 (1995). 



19 



[36] A. Pich, Rep. Prog. Phys. 58, 563 (1995). 

[37] G. Ecker, Prog. Part. Nucl. Phys. 35, 1 (1995). 

[38] R. Koch, Nucl. Phys. A448, 707 (1986). 

[39] L. Alvarez-Ruso, E. Oset and E. Hernandez, Nucl. Phys. A633, 519 (1998). 
[40] G.E. Brown and W. Weise, Phys. Rept. 22, 279 (1975). 
[41] C. Caso et al, Eur. Phys. J. C3, 1 (1998). 



20 



FIGURES 



FIG. 1. Diagrammatic representation of the Bethe-Salpeter equation. 

FIG. 2. Some new terms of the Bethe-Salpeter equation in nuclear matter. Bubbles represent 
particle- hole and A-hole excitations. 

FIG. 3. Im T nn in the scalar isoscalar channel at different nuclear densities as a function of 
the CM energy of the pion pairs. The labels correspond to the nucleons Fermi momentum. 

FIG. 4. Feynman diagrams contributing to the ttN — > ttttN reaction. Dashed lines are pions. 
External solid lines are nucleons. Internal solid lines are nucleons, A's and iV*(1440) where possible. 

FIG. 5. Total cross section for the ttN — > ttttN reaction vs. pion kinetic energy. Upper box: 
tt~p — » 7r + 7r~ra; Experimental points from refs. |26,p7|. Middle box:ir~p — > 7r°7r°n; Data from refs. 



[28H30||. Lower box: ir + p — > 7r + 7r + n; Data from refs. 



FIG. 6. Two pion invariant mass distributions in the it + d —* tt + 7t nn (upper box), and 
7r + + d — > 7r + 7r + (lower box) reactions. Experimental points are from ref. Ji~5|]. 

FIG. 7. Two pion invariant mass distributions in the ir~ + Ca — > it + tt~ X (upper box), and 
7r + + Ca — * ir + 7r + X (lower box) reactions. Solid line, full calculation; dashed line, no medium 
effects in the FSI of the two pions. Experimental points are from ref. ||. 

FIG. 8. Same as fig. [7|, using a fixed averaged density. Short dashed line, p = O.bpo] long 
dashed line p = 0.7po- Solid line, full calculation 

FIG. 9. Same as fig. using a fixed averaged density p = 0.7/?o and a simplified model (see 
text). Dashed line, FSI in vacuum; Solid line, FSI in medium. 
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